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ABSTRACT 

Using a set of multifrequency cross-spectra computed from the three year 
WMAP sky maps, we fit for the umesolved point source contribution. For a 
white noise power spectrum, we find a Q-band amphtude oi A = 0.011 ± 0.001 
/iK^ sr (antenna temperature), significantly smaller than the value of 0.017±0.002 
yuK^ sr used to correct the spectra in the WMAP release. Modifying the point 



source correction in this way largely resolves the discrepancy lEriksen et al.l (120061 ) 
found between the WMAP V- and W-band power spectra. Correcting the co- 
added WMAP spectrum for both the low- l power exces s due to a sub-optimal 



likelihood approximation — also reported by lEriksen et al.l (120061 ) — and the high-£ 
power deficit due to over-subtracted point sources — presented in this letter — we 
find that the net effect in terms of cosmological parameters is a ~ O.Tcr shift in 
ns to larger values: For the combination of WMAP, BOOMERanG and Acbar 
data, we find ns = 0.969 ±0.016, lowering the significance of ng 7^ 1 from ~ 2.7cr 
to ~ 2.0(7. 

Subject headings: cosmic microwave background — cosmology: observations — 
methods: numerical 



1. Introduction 

The results of Wilkinson Mierowave Anisotropy Probe have made an inestimable impact 
on the science of cosmology, highlighted by the very recent release of the three year data: 
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maps, power spectra, and consequent cosmol oKical analysis (jjarosik et al.ll2006l : iPage et al. 
20061 : iHinshaw et al.ll2006l : ISpergel et al.ll2006l ). Precisely because these results play so promi- 
nent a role, it is important to check and recheck their consistency. 



Recently lEriksen et al.l (120061 ) reanalyzed the WMAP three year temperature sky maps, 
and noted two discrepancies in the WMAP power spectrum analysis. On large angular scales 
there is a small power excess in the WMAP spectrum (5-10% at £ < 50), primarily due to 
a problem with the likelihood approximation used by the WMAP team. On small angular 
scales, an unexplained systematic difference between the V- and W-band spectra (few percent 
at £ > 300) was found. In this Letter, we suggest this second discrepancy is at least partially 
due to an excessive point source correction in the WMAP power spectrum. 



Data 



The WMAP temperature data (IHinshaw et al.l |2006| ) are provided as ten sky maps 
observed at five frequencies between 23 and 94 GHz, pixelized using the HEALPijfl scheme 
with 3 million (~ 7'-size) pixels per map. Here we consider the Q-band (41 GHz), V-band 
(61 GHz), and W-band (94 GHz) channels since these have the least galactic foreground 
contamination, but only V- and W-bands for the cosmological parameter analysis. 

We account for the (assumed circularly symmetric) beam profile of each channel inde- 
pendently, adopting the Kp2 sky cut as our mask. This excludes 15.3% of the sky including 
all resolved point sources. To deal with contamination outside the mask, we simply use 
the foreground template corrected maps provided on the LAMBDA websit^. The noise is 
modeled as uncorrelated, non- uniform, and Gaussian with an RMS given by cro,i/ \/Nohs~i{p)- 
Here o"o,i is the noise per observation for channel and Nohs,i{p) is the number of observations 
in pixel p. 



3. Methods 



3.1. Power spectrum estimation 



We estimate power spectra with the pseudo-C^ MASTER method (IHivon et al.ll2002l ). 
which decouples the mode correlations in a noise-corrected raw quadratic estimate of the 
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power spectrum computed on the partial sky. Following iHinshaw et al.l (120031 ) . we include 



only cross-correlations between channels in our power spectrum estimates. 

Considering each of three years, three bands, and the number of differencing assemblies 
per band (two for Q-/V- and four for W-band), 276 individual cross-spectra are available for 
analysis. Each of these is computed to £max = 1024. The V- and W-band spectra have been 
verified against spectra provided by the WMAP team, but the Q-band spectra (computed the 
same way) were not available for comparison. For the point source amplitude analysis, we bin 
the power spectra into ten bins {i = 2-101, 102-201,. . . , 902-1001) in order to increase the 
signal-to-noise ratio and decrease the number of bins (and thus the computation time). The 
corresponding error bars are computed using a Fisher approximation and similarly binned. 



3.2. Point source amplitude estimation 

For our main result, we marginalize over the CMB power and estimate a single ampli- 
tude for the point source spectrum by the method we discuss below. We also compute the 
amplitude in £-bins, but for brevity omit the details, which are similar. We model the ensem- 
ble averaged cross-spectra as the sum of the two components, (Q) = Cl'^^^ + C^'^^'^ , showing 
explicitly the contribution from each part of the signal. Here the multipole bin is denoted by 
/ and the cross-correlation pair by i = {ii){i2) = (Wlyrl)(W2yr3), (Qlyr2)(Vlyr2), etc. No 
auto-power spectra are included, so noise subtraction is unnecessary. We marginalize over 
the CMB spectrum, which we denote by Cp^^. The window functions for each differencing 
assembly pair are w = {wli,}, which we later consider in terms of a matrix. The contribution 
to a cross-spectrum from the CMB signal is thus d''"^^ = w]i,Cp^^ (in thermodynamic 
temperature units). The spectra in this application are already beam- deconvolved, so the 
window functions wji, = are trivial. We denote the amplitude of the unresolved point 
source power spectrum by A. This amplitude relates to the cross-spectra via the frequency 
and shape dependence vector S = {S\}, 

= AS'l 

S\ = w], ir(x(z.,J)(^) K{x{u,,)) 

Ki.) - (1) 
exp(xj 

Here the cross-spectrum i has channels at z/j^ and Ui^, and x(z/) = /izz/fc^TcMB- The units 
of A are antenna temperature squared times solid angle and the function K{x) converts 
from antenna temperature to thermodynamic temperature. Thus, we assume that the radio 
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sources are spatially uncorrelated (and therefore have a white noise spectrum) and have a 
power-law frequency dependence using antenna temperature units. Note that well-resolved 
point sources have already been masked from the maps before the evaluation of the cross- 
spectra, and A therefore represents unresolved sources only. However, we may only directly 
measure the frequency depende nce for the resolved sources. For these, iBennett et al.l (120031 ) 
found P = —2.0, and following iHinshaw et al.l (l2003l . l2006l ) we take the same even for the 
unresolved sources. We choose uq = 40.7 GHz (Q-band) as our reference frequency. 

We organize the binned cross-spectra C/ into a data vector D = {C/}. We use a 
Gaussian model for the likelihood C of the power spectrum, appropriate at high i: 



/:ocexp<'-^[D-(D)]tS-i [D-(D)]) 



(2) 



where the covariance S = ((D — (D))(D — (D))"'') can be written as S = {S",'}. Here we 
assume the covariance is diagonal both in multipole and cross-spectrum. An appendix of 
Huffenberger et al.l (120041 ) der ives an unbiased estirn ator for this type of problem, generalizing 
the point source treatment of iHinshaw et al.l (120031 ). Here the estimators are equivalent, and 
result in a linear estimate for A, denoted A, and its covariance S"^: 



A = (StFS)-iS^FD 
= (S^FS)-\ 

where we have defined the auxiliary matrix 



(3) 



In this notation, we consider D and S as column vectors with a single index il, and w as 
a matrix with indices il and Matrices S and F have indices il and i'l'. This estimator 
marginalizes out the CMB, a conservative treatment which assumes nothing but the fre- 
quency dependence. To compute the amplitude in bins, we redefine A vector of the 
amplitudes, with Cl'^'^'^ = (A • S)], modifying S for each component to lend power only to 
appropriate multipole bins. 



4. Results 

4.1. Point source spectrum amplitude 

Using the method described in the previous section, we find a point source amplitude of 
A = 0.011 ±0.001 /iK^ sr, significantly less than the WMAP value of A = 0.017 ±0.002 //K^ 
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sr (iHinshaw et al.ll2006l ). Computing the spectrum in bins (Figured]), we see that the source 
power spectrum is best measured at 100 < £ < 600. To evaluate goodness-of-fit, we compute 
= J2i hini^i ~ ^y/'^f = 36.6 for 9 degrees of freedom. All of the discrepancy in our fit 
arises from a single high bin at £ = 102-201, which has Ax^ = 27.2. This bin is so different 
that we suspect that it is not detecting point source power alone, but perhaps some residual 
foreground. We leave a rigorous investigation of this anomalous bin to later work, leaving it 
in our analysis here. If we were to exclude it, the other nine bins are consistent with a flat 
power spectrum a.t A = 0.011 /iK^ sr — with = 9-4 for 8 remaining degrees of freedom — 
though they would prefer a somewhat smaller value for A. For the WMAP amplitude, we 
measure = 86.5 for 9 degrees of freedom. This large discrepancy is puzzling because our 
method should be equivalent to the WMAP method. 



4.2. Angular CMB power spectrum 

The net effect of the lower unresolved point source amplitude on the co-added WMAP 
CMB power spectrum may be computed in terms of a weighted average of corrections for 
individual cross-spectra (VxV, VxW, and WxW, respectively). Following the construction 
of WMAP^s spectrum, for i < 500 the correction is given by a uniform average over the 137 
individu al cross-spectrum cor rections; for £ > 500 it is given as an inverse noise weighted 



average (IHinshaw et al.ll2006l ). In this Letter we approximate the latter with the inverse 
variance of the power spectrum coefficients computed from 2500 simulations for each cross- 
spectrum individually, but do not account for correlations between different cross-spectra. 

The net power spectrum correction is shown in Figure [2J We show the Q-band spectra, 
corrected by each point source amplitude, in the top panel of Figure [31 and compare the V- 
and W-band spectra in the bottom panel. 



One of two issues pointed out by lEriksen et al.l (120061 ) was a discrepancy between the 



V- and W-bands at £ > 250 significant at about 3a. This is seen by comparing the two red 
curves in the bottom panel of Figure [31 However, applying the lower point source correction 
raises the V-band spectrum by 10-50 fjK^ in this range but the W-band by only a few /xK^. 
Effectively, about 20 /iK^ of the previous 65 /xK^ average difference is thus removed, reducing 
the significance of the difference from 3 to 2a, compared to 2500 simulations. A small 
difference is still present, and may warrant further investigation, but is no longer striking. 
This gives us confidence that our point source correction is the more consistent than the 
WMAP value. 
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4.3. Cosmological parameters 



To assess the impact of this new high-^ co rrection on cosmological param eters, we repeat 
the a nalysis described by lEriksen et al.l (120061 ) using the CosmoMC package (ILewis &: Bridle 



2002 



which als o gives the paramete r definitions) and a modified version of the WMAP 
likelihood code (IHinshaw et al.ll2006l ). First, at £ < 30 the WMAP likelihood is replaced 



with a Blackwell-Rao G i bbs sampling-ba sed estimator (jjewell et al.l 12004 : IWandelt et al. 



2004 : lEriksen et al.l 12004 : IChu et al.ll2005l ). and second, the bias correction shown in Figure 



[2] is added to the co-added WMAP spectrum. The results from these computations are 
summarized in Tabled! 



As reported by lEriksen et al.l (|2006l ). the most notable effect of the low-^ estimator bias 
in the WMAP data release was a ~ OAa shift in ris to lower values, increasing the nominal 
significance of ris 7^ 1- In Table [1] we see that the over-estimated point source amplitude 
causes a similar effect by lowering the high-£ spectrum too much. Correcting for both of these 
effects, the spectral index is n.s = 0.96 9=fc0.016 for the combination of WMAP, BOOMERanG 



(IMontrov et al.ll2005l : IPiacentini et al.ll2005l : I Jones et al.l 120051 ) and Acbar flKuo et al.l 12004 1 
data, or different from unity by only ~ 2a. The marginalized distributions both with and 
without these corrections are shown in Figure HI The other cosmological parameters change 
little. For reference, the best-fit (as opposed to marginalized) parameters for this case are 
{%h^, rich"^, h, T, ns, MlO^Mg)} = (0.0225, 0.108, 0.732, 0.919, 0.967, 3.05). 



5. Conclusions 

Using a combination of cross-spectra of maps from the Q-, V-, and W-bands of WMAP 
three year data, we fit for the amplitude of the power spectrum of unresolved point sources 
in Q-band, finding A = 0.011 ± 0.001 /iK^ sr. This fit has significantly less power than the 
fit used to correct the WMAP final co-added power spectrum used for cosmological analysis. 

We compute and apply the proper point source correction, noting the corrected V- and 
W-bands are more consistent than before. The improper point source correction conspires 
with a low-i estimator bias to impart a spurious tilt to the WMAP temperature power 
spectrum. With the revised corrections, we find evidence for spectral index ns 7^ 1 at only 
~ 2cr, while other parameters remain largely unchanged. 
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Table 1. Cosmological parameters 



Parameter WMAP Low-£ and PS corrected 



WMAP data only 





0.0222 ± 0.0007 


0.0223 ± 0.0007 




0.241 ± 0.036 


0.244 ± 0.035 


log(lOiMs) 


3.019 ± 0.067 


3.039 ± 0.068 


h 


0.731 ± 0.033 


0.730 ± 0.032 


ns 


0.954 ± 0.016 


0.966 ± 0.016 


T 


0.090 ± 0.030 


0.090 ± 0.030 


WMAP + Acbar + BOOMERanG 




0.0225 ± 0.0007 


0.0225 ± 0.0007 




0.239 ± 0.031 


0.240 ± 0.031 


log(lOiMs) 


3.030 ± 0.064 


3.045 ± 0.065 


h 


0.737 ± 0.029 


0.738 ± 0.030 


ns 


0.958 ±0.016 


0.969 ± 0.016 


T 


0.091 ± 0.030 


0.091 ± 0.030 



Note. — Comparison of marginalized parameter results 
obtained from the WMAP likelihood (second column) and 
from the WMAP + Blackwell-Rao hybrid, applying the low- 
i estimator and high-£ point source corrections (third col- 
umn). 
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Fig. 1. — The point source power spectrum, fit from WMAP Q-, V-, and W- bands. Tlie 
lowest i bin is not plotted, because tfie error bars span tlie entire range of the plot, and it 
has little statistical influence. The 1-par ameter fit for a flat sp ectral shape is also shown, as 
well as the point source amplitude from iHinshaw et al.l ( 20061 ). 
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Multipole moment, / 

Fig. 2.— The net difference [{C^""^ -Cf'^)^{^+l) /2ti] in tfie final co-added WMAP spectrum 
due to the new and smaller point source amplitude. The sharp break at £ = 500 is due to 
different weighting schemes, and the smaller fluctuations at high £'s are due to a finite number 
of Monte Carlo simulations for noise estimation. 
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Fig. 3. — The impact on the Q-, V-, and W-band power spectra of the revised point source 
correction. The top panel shows the Q-band power spectrum with the WMAP point source 
correction (red) and the correction in this work (green dashed), plotted with the WMAP 
best-fit ACDM spectrum. Particularly at i < 400 where noise is lower, this highlights the 
point source over-subtraction using the WMAP correction. The bottom panel shows the V- 
band (soli d) and W-band (dash ed) power spectra minus the co-added WMAP temperature 
spectrum (IHinshaw et al.ll2006l ). computed with the WMAP point source correction (red) 
and the correction in this work (green). The V- and W- bands are internally more consistent 
with the revised source correction. 
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Fig. 4. — Marginalized posterior distributions for the spectral index ng computed with the 
combination of WMAP, BOOMERanG and Acbar data, both with the WMAP likelihood 
code as provided (dashed) and after applying a low-£ estimator correction and a high-£ point 
source correction (solid). 



